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Abstract 

The CUJET Monte Carlo Jet Energy loss model is applied to predict the jet flavor, centrality 
and density dependence of the nuclear modification factor Raa at RHIC and LHC. Running 
coupling eff'ects due to combined x, and evolution are included for the first time in the 
dynamical DGLV opacity expansion framework and are shown to provide a natural dynamical 
QCD tomographic solution to the surprising transparency of the quark gluon plasma produced at 
LHC as suggested by > 10 GeV Raa data from ALICE, ATLAS, and CMS. 



1. Introduction 

In this proceeding we report results of a recent update of the CUJET d pQCD tomographic 
model, where running coupling eff'ects have been taken into account. We remind that CUJET 
is a model developed as part of the ongoing Department of Energy JET Topical Collaboration 
121 effort to construct more powerful numerical codes, which are necessary to reduce the large 
theoretical and numerical systematic uncertainties that have hindered so far quantitative jet to- 
mography. CUJET extends the development of the the GLV, DGLV and WHDG 1 3 , 4 , 5 1 opacity 
series approaches by including several dynamical features that allow the evaluation systematic 
theoretical uncertainties such as sensitivity to formation and decoupling phases of the QGP evo- 
lution, local running coupling and screening scale variations, and other eff'ects out of reach with 
analytic approximations. 

The model is here applied to predict the nuclear modification factor RAA^a^fiy ~ 0, , ^/s, C 
- 5%), for a variety of jet parton flavors a = g,u,c,b and final fragments f = 7i,D, B, over 
a broad pr kinematic range at mid-rapidity for central collisions at = 0.2, 2.76 GeV as ob- 
served at RHIC |6| and LHC |7|. By implementing running coupling eff'ects, justified by the 
broad range of energies at play and already suggested in 1 8 1, we want to show that this relatively 
simple physical phenomenon could account for the surprising transparency of the QGP at LHC, 
as recently observed by the ALICE and CMS experiments |i9l . 

2. CUJET model 

The Monte Carlo techniques implemented in CUJET allow the model to compute finite order 
in opacity n > I contributions to the jet medium induced gluon radiative spectrum. As in previous 
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works fPl, however, we limit our studies to the first order in opacity; furthermore, we model the 
interaction potential between the jet and the medium with the pure HTL dynamical model |10|. 
Neither of these approximations are shown to qualitatively alter the results, once rescaling of the 
eff'ective coupling constant — or as we will see shortly — is taken into account. Both radiative 
and elastic contributions to the energy loss are considered, and fluctuations of the radiated gluon 
number included via Poisson expansion (incoherent emission). 

The plasma is assumed inhomogeneous (Glauber profile) and nonstatic (1 + lD Bjorken ex- 
pansion), and its density profile is constrained solely by the initial observed rapidity density 
dN/dy, equal here to 1000 for = 0.2 GeV RHIC colHsions and 2200 for = 2.76 GeV 
LHC collisions; we also set the initial thermalization time tq = 1 fm/c and take T 100 MeV as 
the characteristic temperature at which the jet decouples from the medium. 

3. Running coupling 

As we already mentioned, recent preliminary results published by the ALICE and CMS col- 
laborations | 9 | appear to indicate a steeper rise of Raa in the range of momenta 30-100 GeV. In 
addition to this, LHC extrapolations performed using an eff'ective static coupling constant fit to 
RHIC results seem to systematically overpredict the quenching suff'ered by the jets in the plasma 
ifm . These findings motivated us to relax the eff'ective fixed alpha approximation and utilize a 
one-loop order running coupling, parametrized as follows 1 12 |: 



Here we choose to keep ao as the only free parameter of the model. The choice of scale Q, of the 
order of 1 GeV, is somewhat arbitrary: of the three powers of as that occur in the DGLV integral, 
two are related to the medium scattering vertex and one to the radiated gluon. We let the first two 
to scale as q^, the transverse momentum exchanged between the jet and the medium, and the 
latter to scale approximately with t ^ k]_/(x(l- x)), where k represent the transverse momentum 
of the radiated gluon and x is its fraction of plus-momentum. In order to account for this source 
of systematic uncertainty, we let the scale vary of ±25%, while fixing the parameter to fit one 
chosen pion R^T^Pt = ^OGeV) = 0.35 point. 

At the same time, we also include running coupling eff'ects in the elastic contribution to the 
energy loss 1 13 |. The two powers of as, in this case, scale as ET and yu^, with E,T,/d the energy 
of the jet, the temperature of the plasma and the Debye screening mass respectively. 

Of all these contributions, the most relevant comes from the radiated gluon vertex, as(\^]_ /(x(l- 
x))); in particular, this helps explaining the faster rise in the shape of pion Raa- In fact, the hard 
tails of the opacity expansion spectrum make the integral, whose upper kinematic limits are a 
function of the energy of the jet E, more sensitive to the high k^ region where as is smaller. This, 
in turn, clearly modifies the energy dependence of the nuclear modification factor by reducing 
the relative energy loss at large energies; as a consequence, Raa will show a steeper rise with 
than in the fixed coupling scenario. 



In CUJET, the parameter ao in the expression for the running coupling is the only free pa- 
rameter that we fit to the data: constraining one reference pr = 40 GeV point of pion R^^ = 0.35 
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4. Results 




Figure 1: (Left) Pion Raa at LHC (solid lines) and RHIC (dashed lines). In black the fixed coupling CUJET results, 
constrained at pr = 10 GeV RHIC with as = 0.3. In red the running coupling CUJET results, constrained at pr = 40 
GeV LHC with ao = 0.4. Central 0% - 5% preliminary ALICE and CMS h± LHC data 1 9 1 (brown and gray triangles, 
respectively) are compared to predictions. (Right) Notice the accordance with data [I4j of the RHIC Pion Raa running 
coupling result in the range of energies 5-20 GeV. 



at LHC, we set ao = 0.4. The backward extrapolation to RHIC is then parameter free, assuming 
that ao does not vary with ^/s. This is an inversion of the standard practice to fit the data at 
RHIC and extrapolate to LHC, however is needed since the physics we now want to probe spans 
over a much broader range of energies, inaccessible at RHIC. The results are shown in Fig. 1. 
Observing the figure on the left, it is evident that the overall shape of Raa across the broad range 
of Pt under consideration is changed with respect to the previous fixed coupling results. Besides 
appreciating the more satisfactory agreement with the data, both at LHC and RHIC (in the latter 
case our predictions are almost left unchanged given the restricted range of energies at play, as 
shown in the panel on the right), it is surprising to note how the eff'ective energy dependence 
itself of the energy loss appears to be modified. 

In Fig. 2, we are showing the flavor dependence of Raa and in particular the comparison 
with the preliminary D meson data from ALICE ||9l . It is worth restating that the heavy meson 
predictions plotted on the left are constrained by the same parameters already fit to the pion 
results, a fact that seems to indicate a remarkable robustness of CUJET, at least when compared to 
the presently available data. On the right panel, we can instead observe the full flavor dependence 
of the nuclear modification factor: the most striking feature is the inversion of the n < D < B 
Raa hierarchy ordering at high pr, already seen in the case of the fixed coupling constant in HI. 
This eflfect is mostly due to the steeper initial invariant jet distributions of c and b jets at RHIC 

na. 



5. Conclusions 

The CUJET model has been applied here to study the flavor and ^[s dependence of the 
nuclear modification factors for central collisions at mid-rapidity. The inclusion of running cou- 
pling eff'ects in the model has vastly improved the agreement with the recent LHC data without 
aff'ecting the RHIC sector in the range of energies currently probed. With one free parameter 
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Figure 2: (Left) D meson Raa at LHC. In black the running coupling CUJET results, constrained by the Pion fit shown 
in Fig. 1. Central 0% - 7.5% preliminary ALICE LHC data | 9 | are compared to predictions. (Right) Illustration of jet 
flavor tomography level crossing pattern of central Raa versus pr for Pions, D and B mesons. LHC Pb+Pb predictions 
are shown in solid color, RHIC Au+Au results are shown in faded colors. The opacity is constrained at LHC, given 
dN/dy(LHC) = 2200, by a fit to a reference point R\^au^Pt = 40 GeV) = 0.35 setting ao = 0.4. 



(ao) used to fit the pion data at LHC, we have predicted the same level crossing pattern of Raa 
already seen in the case of fixed coupling. At the same time, we have shown that the theoret- 
ical uncertainties relative to the choice of running scales do not aff'ect qualitatively the results. 
Further data, especially for heavy mesons (both at RHIC and LHC) and extended energy range 
(RHIC) may give additional indication of running eff'ects and lead to stronger conclusions with 
respect to the energy dependence of the energy loss itself. 
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